Quantum-circuit guide to optical and atomic interferometry 
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Atomic (qubit) and optical or microwave (modal) phase-estimation protocols are placed 
on the same footing in terms of quantum-circuit diagrams. Circuit equivalences are used 
to demonstrate the equivalence of protocols that achieve the Heisenberg limit by employing 
entangled superpositions of Fock states, such as NOON states. The key equivalences are 
those that disentangle a circuit so that phase information is written exclusively on a mode 
or modes or on a qubit. The Fock-state-superposition phase-estimation circuits are converted 
to use entangled coherent-state superpositions; the resulting protocols are more amenable to 
' realization in the lab, particularly in a qubit/cavity setting at microwave frequencies. 
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Atomic and optical (or microwave) interferometers are standard tools for high-precision metrol- 
ogy: the quantity to be measured is mapped onto the phases of internal levels of atoms or the 
paths or polarizations of photons, and the phase is estimated by interferometry. The conventional 
C<") | way of operating an interferometer, using independent (unentangled) atoms or independent (un- 

entangled) photons, i.e., photons created and processed by the techniques of linear optics alone, 
has a sensitivity limit that scales as l/y/N, where N is the number of atoms or photons devoted 
to the task. This scaling is called the quantum noise limit or the shot-noise limit 

The sensitivity of an interferometer can be improved by entangling the atoms or photons. The 
ultimate sensitivity is obtained by using so-called cat states in atomic interferometry or 
their analogue, NOON states, in optical interferometry [j| @, Q, @, 0, 0, 11, G3]. The moniker 
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NOON, which comes from the (unnormalized) form, \N, 0) + |0,iV), was applied and popularized 
by Dowling and collaborators jl, Q3]. The use of entanglement permits a sensitivity that scales 
as 1/N, a scaling often called the Heisenberg limit [lil . 14 . [l5| . There are many closely related 
ideas for achieving Heisenb erg-limited scaling, some of which don't look like an interferometer at 
all (because they aren't) and don't require any entanglement. In this paper we use the pictorial 
representation of quantum-circuit diagrams to put atomic and optical interferometry on the same 
footing and to demonstrate the equivalence of various methods for achieving the Heisenberg limit. 

In its translations from atoms to photons, this paper formulates a quantum-circuit version of 
the equivalence between atomic and optical interferometry demonstrated by Yurke, McCall, and 
Klauder [la ]. Much of the paper's content onphase estimation and interferometry is anticipated 
by the work of Gerry and collaborators [1, 0, Ha, 0, [13] , and many of the qubit-mode operations 



discussed in the paper follow a path blazed by Davidovich, Haroche, and collaborators [la, LL9|, l2d . 
ELS [23], [24|. Lee, Kok, and Dowling [§| made an initial foray into the domain of quantum circuits 
as tools for investigating equivalent atomic and optical interferometers, calling this the "quantum 
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Rosetta stone." Less lofty, but more ambitious, this article aims to establish a permanent beachhead 
on this terrain, making it a standard tool for analyzing interferometers and phase estimation. 

This article is dedicated to the memory of Krzysztof Wodkiewicz. When new ideas came along, 
as in the quantum-information revolution, Krzysztof embraced them and incorporated them into 
his suite of theoretical tools. Indeed, this contribution is very much in the spirit of his own 
recent work on qubit decoherence [25], [2^, 27, 0, [2^], something he knew as well as anyone, but 
which he rethought, reexamined, and reformulated in the light of the modern language of quantum 
operations and Kraus operators. In the same way, this contribution repackages concepts from 
optical and atomic interferometry, concepts that all physicists know, but presented here from a 
different perspective. It is to be hoped that the paper lives up to Krzysztof 's high standards. The 
paper also represents an experiment that CMC has long wanted to try: writing a paper in which 
the figures and captions outweigh the text. 

Section |Tl] describes atomic interferometry in terms of qubits and optical interferometry in terms 
of modes, reviewing how to map the states and unitary operations for qubits to those for modes; 
it introduces the language of quantum circuits, with particular attention to how measurements 
are represented and manipulated within the quantum-circuit picture and the role of disentangle- 
ment in interferometry. Section IIIII reviews the description of conventional, quantum-noise- limited 
interferometers in terms of circuit diagrams. Section HVl turns to Heisenb erg-limited modal inter- 
ferometers. It considers phase-estimation protocols that are best thought of as using superpositions 
of Fock states, and Sec. |V] discusses the analogous protocols that use superpositions of coherent 
states. A concluding Sec. IVII briefly considers how Heisenb erg-limited phase estimation might be 
implemented in the lab and then wraps up with a peroration to circuit diagrams as the best way 
to represent general protocols for phase estimation and interferometry. 



II. QUBITS AND MODES 
A. States, gates, and quantum circuits 

Atomic interferometry deals with atoms, idealized as having two levels. Any such two-level 
system can be regarded as a qubit, which has two standard states, a "ground" state |1) and an 
"excited" state |0). Thus we often refer to atomic interferometry as qubit interferometry. 

The pure states of a qubit are conveniently thought of as lying on a Bloch sphere defined by 
Pauli operators Z = a z = |0)(0| - |1)(1|, X = a x = |0)(1| + |1)(0|, and Y = a y = -i(|0)(l| - |1)(0|). 
We use Latin letters at the beginning of the alphabet to represent bit values, and 1, and we 
use Latin letters at the end of the alphabet to denote the sign representation of a bit value, 
e.g., Z\a) = (— l) a |a) = z\a). When we have N qubits, we introduce a total "angular momentum" 
operator J = \ X^=i a r We denote the product basis of standard states by \a) = \a\) <g> ■ ■ ■ (g> |ajv), 
where a stands for the bit string a\ . . .a^. For qubit interferometry, we are interested in symmetric 
qubit states. The (N + l)-dimensional symmetric subspace is spanned by the states 



no, ni) = I J = N/2, m = (no — n\)/2) = 




a 



where no is the number of atoms in state |0) and n\ is the number in state |1), and where the sum 
runs over all strings a that have no Os and n\ Is. 

In optical or microwave interferometry we deal with two modes of the electromagnetic field, 
labeled and 1, with annihilation operators a and b. Thus we often refer to optical interferometry 
as modal interferometry. The number operators for the two modes are denoted Ao = a)a and 
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JVi = tfb. The Fock basis of photon-number eigenstates for the two modes is defined by 



(a 



no,wi)~ = K ' K ' \0,0)~ , (2.2) 



where |0, 0)^ is the vacuum state of the two modes, and where no is the number of photons in mode 
and n\ the number in mode 1. The connection to qubit interferometry is made by identifying these 
Fock states with the symmetric qubit states denoted in the same way. The squiggle in Eq. (|2.2[) is 
the symbol we use for modes; we almost always omit it in state designations. For N = 1, the Fock 



states are the dual-rail logical states of linear-optical quantum computing [30|, l3l(, |1,0)^ = |0) 
and |0,1)~ = |1). For N = 2, the Fock states are |2,0)~ = |00), |1,1)„ = (|01) + \W))/V2, and 
|0,2}„ = |11). 

The identification of symmetric qubits with modes is completed by introducing the angular- 
momentum operators of the Schwinger representation: 



J z = ^(cJa-tfb) = ^J2z j , (2.3a) 

3=1 

1 1 - 

4=2^ + ^) = -^;, (2.3b) 

3=1 
1 N 

J y = -U<Jb-tfa) = -Y J Y j . (2.3c) 



3=1 

The photon-number operator is N = a) a + b^b. The identification of the iV-photon subspace of 
the two modes with the symmetric subspace of N qubits is the foundation for unifying optical and 
atomic interferometry. Of course, when dealing with modes, we can also have superpositions of 
different total photon numbers. These superposition states do not have a convenient representation 
in terms of qubits, but we will be able to represent them in our circuit diagrams. 

The allowed unitary operations, or gates, are symmetric unitary operators on the N qubits or, 
equivalently, photon-number preserving unitaries for the two modes. All such unitaries, up to a 
global TV-dependent phase, can be written as e~ i ^ J \ where /(«/) is a Hermitian function of the 
total angular momentum [321 ] . 

An important special case occurs when / is linear, in which case the unitaries, rotations 
e -ij-nB _ ^ e -io--n9/2^®N j n an g u l ar . m0 mentum sense, are independent, identical unitaries 
on the qubits or, equivalently, the operations of linear optics on the two modes. The operation of 
such a unitary is summarized by how it transforms the modal creation operators: 

e -*j.<w ( aj \ e% j.™ = m t ( aj \ (2 4) 



Here M is a 2 x 2 unitary matrix independent of N, and T denotes the transpose. Applying both 
sides of Eq. (|2.4p to the modal vacuum state and converting to qubit notation gives 



e 



from which it follows that M is the representation of e l(T " ne / 2 in the standard qubit basis, i.e., 
M a b = {a\e~ teT " n6 / 2 \b). The matrix M is given explicitly by 

= / cos(d/2) - in z sm(0/2) {-in x - n y ) sin(6*/2) \ . . 

\ (-in x + n y )sm(0/2) cos(^/2) + in z sin(fl/2) J ' [ 1 
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FIG. 1: Linear-optical quantum gate in qubit form on the left, in terms of Schwinger operators for qubits 
in the middle, and in modal form on the right. The qubits provide a particle description, and the modes 
provide a wave description: the qubits are the particles — say, atoms or photons — that occupy the modes 
and 1, and the modes are the single-particle states, |0) and |1), of the qubits. In qubit circuits we use 
N = 3 qubits for illustration (occasionally N = 5), here and in the following, but the equivalences to modal 
circuits are valid for arbitrary N. A state to the left of a circuit wire indicates the input to that wire; if no 
state is given, equivalences are independent of the input state, as they would be here. Equivalences between 
qubit and modal circuits require, however, that the input state be in the symmetric subspace of iV qubits 
or, equivalently, in the TV-photon sector for modes. On the left and right, wires are labeled explicitly by 
whether they are qubits or modes; hereafter, we omit this labeling and use a squiggle, as in the right circuit, 
to designate a modal wire. When n lies in the equatorial plane, these are, for 8 = tt/2, tt/2 transitions for 
qubits and 50/50 beamsplitters for modes and, for 9 — n, qubit FLIP operators and modal SWAP operators. 



Suppose n = x cos (ft + y sin</> lies in the equatorial plane. If 6 = tt/2, these are ir/2 transitions 
for qubits and 50/50 beamsplitters for modes, with transformation matrix (|2,6p given by 

M=-^=( A* "TH • (2-7) 



y/2 V 1 

If 9 = 7r, these are qubit FLIP gates and modal SWAP gates, with additional phase shifts; i.e., 
e -iJ-n-K _ . fi)® N has transformation matrix 

m =(-«•* if*)- < 2 - s > 

Notice also that Z® N = i N e~ lJz7T = e jfttb?r = / <g> II, with LT being the modal parity operator, and 
that (-Z)® N = (-i) N e-' iJz7T = e - iata7T = 11®/. A 2vr rotation about any axis h has M = I and 

e -u-h2n = (_-giv = n®n. 

These considerations provide an opportunity to introduce quantum circuits in Figs. [Q and El 



We use use generally conventional circuit notation [33|], with wires representing qubits or modes 
and unitary gates enclosed in rectangular boxes, but we use an initial squiggle to indicate when a 
wire applies to a mode. We make frequent use of the qubit Hadamard and S gates: 



H = ie -i[{x+Z)/^/2 = _L {X + Z) ^ _L f J _M 



(2.9) 



S = e iir/4 e -tZir/4 ^ (j i) ' (2 ' 10) 

Notice that e~ lY ^^ = (I—iY)/y/2 = XH = HZ. Controlled application of a unitary is represented 
conventionally by a circle on the control qubit connected by a line to the unitary on the target. 
Control on |1) is represented by a filled circle, and control on |0) by an open circle. Figure [1] is the 
circuit form of the equivalence between qubit and modal gates for linear-optical transformations. 
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FIG. 2: Optical Kerr gates and their qubit (particle) equivalents: (a) a controlled-SIGN (controlled-Z) gate 
between every pair of qubits is equivalent to an ordinary tt/2 phase shift on mode 1 plus a tt/2 self-Kerr 
phase shift on mode 1, i.e., a phase shift due to a coupling proportional to the square of the number of 
photons; (b) adding Z rotations and an overall phase shift to (a) gives a cross-Kerr coupling between modes; 
(c) a Z-rotation by angle N9 on each qubit is equivalent to Kerr phase shifts on each mode. As noted in the 
text, these gates arise naturally in discussions of nonlinear interferometry based on two-particle (quadratic) 
interactions; circuit (a) appears again in Fig. I2"27c) as part of a circuit for making mode-entangled coherent 
states. The equivalence in (a), which is essentially obvious, can be regarded as coming from the Hamiltonian 
identity H = \ E J »fe(^- / )(^- / ) = \Jt~\{N ~l)Jz + \{N 2 /2-N) = ±Ni(iVi-l); e~ iHv is a controlled- 
SIGN gate between every pair of qubits. The equivalence in (b) combines this Hamiltonian identity with 
^NqN\ — \N 2 — |J|; (c) comes from \{Nq — N 2 ) — NJ Z . Since the gates in these equivalences are all 
diagonal in the standard basis, the equivalences can be mixed without regard to ordering of the gates. 




FIG. 3: Photon counting on a mode, i.e., a measurement of the photon-number operator a^a, is represented 
by the conventional symbol for a photodetector. 



Figure [2] gives examples of gates and equivalences outside the linear-optical set; these equivalences 
lie at the heart of recent work on nonlinear modal interferometry Hi, EE H EH 13. 

The qubit and modal circuits of Figs. [T] and [5] illustrate strikingly the difference between en- 
tanglement among qubits, which we call particle entanglement, and entanglement between modes, 
which we call wave entanglement. The qubits at the output of the circuit in Fig. [T] are obviously 
not entangled, yet in the case of a 50/50 beamsplitter, for example, the output modes are entan- 
gled: the two modes have anticorrelated, binomial photon-number distributions in superposition. 
A beamsplitter is an entangling gate for modes, but it is a product gate for qubits. In Fig. [21 
the entangling controlled-SIGN gates in (a) are equivalent to a product gate for modes; the gates 
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in (b), which for qubits differ from those in (a) only by a set of local gates, entangle both qubits 
and modes; the gates in (c) are products for qubits and for modes. 

We emphasize that qubit and two-mode interferometers can always be thought of in both ways, 
as involving particles or modes — that's the complementarity of quantum mechanics — but we prefer 
to think one way or the other, perhaps to our detriment. We generally like to regard atomic 
interferometry as using particles and optical or microwave interferometry as using modes, but the 
other way of thinking is available in both cases. There is one situation, atomic interferometry using 
a two-component Bose-Einstein condensate (BEC), where it is so natural to think in both ways 
that it is especially valuable to be able to translate between the two ways of thinking. 



B. Measurements and circuit equivalences 

The remaining ingredient needed in our circuit diagrams is a representation of measurements. 
For qubits, we are interested in symmetric measurements and thus can specialize to measurements 
of a particular Pauli component a • n, usually Z, on all qubits. We write the ±1 outcome of a 
measurement of a Pauli component as z = (— l) a . Most discussions of qubit measurements regard 
the bit value a as the outcome, but in this paper we always think of z as the outcome. For modes, 
the only measurement we need is the measurement that counts photons in a mode. 

In our circuit diagrams, measurements are represented by bubbles (rounded boxes), with the 
measured observable inside the bubble. For photon counting, we adopt a special notation, il- 
lustrated in Fig. [3l in which the measurement is represented by the conventional symbol for a 
photodetector. We are only interested in the outcome of a measurement and its statistics and not 
in the post-measurement state of the measured system; in the lingo of quantum measurement the- 
ory, this means we are interested in the measurement's POVM, but not in the quantum operation 
that outputs a post-measurement state. Thus, in our circuits, there is no quantum wire emerging 
from a measurement bubble; the only output from a bubble is a double (classical) wire that carries 
the value of the outcome. The values on such classical wires can be manipulated by classical gates, 
which we depict as rectangular boxes. We often label a classical wire with the value it carries. 

Figure H] illustrates our measurement notation in the context of two ways of processing the data 
from Z measurements on all qubits and the equivalent measurements on modes. The important 
conclusions are that summing the outcomes of Z measurements is equivalent to counting pho- 
tons on the two modes and then differencing the counts, whereas multiplying the outcomes of Z 
measurements, i.e., measuring Z® N , is equivalent to measuring the parity of mode 1 aa a. 

We need several additional rules for handling measurements, the first three of which, nota- 
tion for a measurement whose outcome is discarded, the rule for unitarily transforming measured 



observables, and the principle of deferred measurement [33|, [37(, are introduced and described in 
Fig. [5J Figure [6] combines the equivalences of Fig. [5] to derive a new set of equivalences, spelled out 
in Fig. [6jb) , which lie at the heart of interferometry. These equivalences express a very important 
concept, disentanglement, in several different ways. If there is a most important figure to master 
before considering Heisenb erg-limited interferometers, Fig. [6] is it. In using the equivalences of 
Fig. EJb), it is important to remember that they require that the unitary operator U be Hermitian, 
which means that it is an observable all of whose eigenvalues are equal to ±1, and that for the last 
two circuits, the ±1 eigenspaces of U have the same dimension. 

The equivalence of Fig. [6ja) is also used in Fig. Etc) to derive a standard method for measuring 
a Hermitian U, of which the parity of a field mode is an outstanding example. This measurement 
is founded on a technique called phase kickback [H, which maps the phase of U, ±1 in this 
case, to a qubit. Phase kickback is a key element in many quantum algorithms. This qubit-assisted 



method for measuring the parity of a mode was proposed by Lutterbach and Davidovich [2l|, [2 
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FIG. 4: (a) Measurement of Z on TV = 3 qubits, all of which are prepared in the state (|0) + \l))/V2, 
followed by classical processing that outputs the sum of the TV outcomes. The final outcome 2m = no — n± 
is the difference between the numbers of qubits found in |0) and |1); it is distributed binomially, with mean 
zero and variance TV. The second and third circuits transform from qubits to modes via the Schwinger 
operators, with the third circuit having TV photons in mode processed through a 50/50 beamsplitter and 
then subjected to a measurement of the photon-difference operator 2J Z — Nq — N\. In the final circuit, 
the differenced photocounting is achieving by photocounting on each mode, followed by differencing the 
counts to give outcome 2m = Uq — n\. (b) Measurement of Z® N on TV = 3 qubits prepared as in (a), with 
outcome (— l)™ 1 equally likely to be ±1, is equivalent to measuring Z on all qubits, followed by classical 
processing that outputs the product of the TV outcomes. This measurement is equivalent to measuring the 
parity II = (— l) b b of mode 1. In the final circuit, the parity measurement is achieved by photocounting 
followed by classical extraction of the parity. Since mode is unmeasured, it can be restored as a quantum 
wire in the final circuit even though this is not equivalent, strictly speaking, to the third circuit. The quoted 
probabilities for the outcomes depend on the initial state shown, but the circuit equivalences themselves are 
independent of initial state. 



and implemented in Refs. [39|, [40|. 

The equivalences in Fig. [6^b), though easy to derive from circuit equivalences, have wandered 
far enough down the road of circuit diagrams that it is instructive to check that we have not gone 
astray. A good way to do this is to demonstrate the equivalences in the conventional language of 
state vectors. Suppose that in Fig. Mjo), the input state has the relative-state decomposition 



\f) = vWo> 8) |0) + VWi) ® |1) , 



(2.11) 
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FIG. 5: (a) Measurement of an observable O, whose outcome is discarded, is represented by a bubble without 
a classical wire emerging from it. Discarding the outcome implies averaging over the possible outcomes; this is 
equivalent to tracing out the measured system. We can always represent this as a no-outcome measurement, 
i.e., a measurement of the unit operator /; since a measurement of / is really no measurement at all — it merely 
instructs us to trace out that system — we can choose to restore the system to the circuit by eliminating the 
measurement of I and extending the quantum wire, provided we remember that the wire must ultimately be 
terminated by a measurement without interacting with other wires in the circuit, (b) Top equivalence: Rule 
for unitarily transforming a measured observable. Preceding a measurement of observable O — Xk\k)(k\ 
with unitary U is equivalent to measuring the conjugated observable WOU = Ylk ^kU^\k)(k\U; in both 
cases, incoming state CkU'\k) yields outcome with probability \ck\ 2 ■ Bottom equivalence: Special case 
in which P is a unitary that permutes the eigenstates of O, i.e., P\k) — \P(k)), where P(k) is a permutation 
of the indices. In this situation, P^OP — J2k ■^k\P~ 1 (k)) (P~ 1 (k)\ = Xp(k)\k){k\; thus a measurement 
of P^OP is the same as measuring O and then applying the permutation to the outcome. On both sides of 
the equivalence, incoming state J2k c k\k) yields outcome Xp(k) with probability \ck\ 2 ■ (c) Top equivalence: 
Principle of deferred measurement, depicted here for a mode and a qubit. Measuring Z on the qubit and 
using the outcome to control classically the application of a unitary U to the mode is the same as applying 
a controlled-?/, with the qubit as control and the mode as target, and then measuring Z on the qubit. In 
the latter case, the control is done coherently, and the outcome of the measurement later reveals whether 
the unitary was applied. Bottom equivalence: Same principle, but with the outcome of the Z measurement 
discarded after application of the classical control. 

where the states |V>o) and \tpi), though normalized, are not necessarily orthogonal, and qo + qi = 1. 
We show that circuits (i)-(iv) of Fig. E^b) all lead to the same probability for the outcome z. 
The controlled-?/ in circuit (i) transforms the input state to 

VWo>® |0} + v^Wi>®|1> , (2.12) 
and the subsequent Hadamard gate, H, leaves the state 

-4(vWto> + V9i£W) ® |0} + 4=(x/9o|^o)-VgI?7|Vi}) , (2-13) 





FIG. 6: The unitary operator U is also Hermitian, i.e., U = U', implying that all its eigenvalues are equal 
to ±1. (a) Conjugation of / ® X by a controlled-{7 gives U ® X . In the middle circuit, X has been pushed 
through the initial controlled-?/, turning the control on |1) to control on |0). (b) Circuit equivalences for 
disentangling circuits. The first two steps unitarily transform the measurement of /<g> Z to U ® X. The next 
step reverts to a Z measurement preceded by a Hadamard and includes explicit classical control circuitry 
for multiplying the outcomes of the U and Z measurements. The fourth step uses the bottom equivalence 
in Fig.[5^b) to replace the post-measurement classical multiplication, which swaps ±1 outcomes, with a pre- 
measurement unitary P that swaps the ±1 eigenspaces of U, i.e., PUP^ = — U . This step thus requires that 
the ±1 eigenspaces have the same dimensions. The final step uses the principle of deferred measurement of 
Fig. [^c) to move the controlled operation through the Z measurement. The labels (i)-(iv) are introduced 
for later reference. That (i)-(iv) all lead to the same probability for the ultimate outcome z is worked out in 
detail in the text, (c) Measurement of U, achieved with the use of an ancillary qubit. The first equivalence 
is true because the qubit's input state is |0). The second equivalence follows from using Z — HXH twice, 
surrounding an application of the equivalence in (a). The result is a standard method for measuring parity II. 
Notice that this part does not require that U have ±1 eigenspaces of the same dimension. 
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from which it follows that the probability to obtain outcome z in the Z measurement on the qubit is 



Pz 



golV'o) + Zy/qiU\ipi 



2 = \{ l +2: v ^o?TR«>((r„|rr|r l }) 



(2.14) 



For circuit (ii), the probability for outcome z follows from the expectation value of U <S> X in 
the input state (|2.1ip . i.e., 



Pz 



(2.15) 



which is the same as Eq. (|2.14p because U = U* . 

In circuit (iii), the state after the initial Hadamard takes the form 



v/2 



= (vWo> + vWi>) ® 1°) + >/^;^=(VWo) " vWi>) ® |1) , (2.16) 



where 



golV'o) +yy/qi\ q Ph 



2 1 



l + 2y^/^TRe((V'o|^i)) 



(2.17) 



is the probability for the Z measurement to yield outcome y. The state of the mode after the Z 
measurement yields outcome y is 



1 



(2.18) 



Hence, the conditional probability for the measurement of U to have outcome x, given outcome y 
for the Z measurement, is 

1 / / expectation value\ \ 



'y - •y 
The unconditioned probability for the outcome of the U measurement is 

Px = z^p x \y r y = ^(1 + ^(90(^01^1^0) + qii^iW^i, 



(2.19) 



(2.20) 



In contrast, the conditional probability for the product of the two outcomes, z = xy, is p z \ y 
Px=zy\yi which gives an unconditioned probability 



Pz = ^2Px=zy\ y r y = 2 (l + 2z^/q^Re((ip \U\^i))^ , 



(2.21) 



the same as that in Eq. (|2. 141) . 

In circuit (iv), after the initial Hadamard and the controlled-P, the state becomes 

-4(v«biVto>+vWi>) ® \^)+\ p (yTo\^)-v^\^)) ®ii> • (2-22) 



11 



Hence, the probability of outcome z in the measurement of U is 

= ~ (l + 2z^mM^o\U\M) , (2.23) 

where the final step uses that P conjugates U to its opposite, i.e., PUP^ = —U. The final 
probability is again the same as that in Eq. (|2.14p . 

There is an important special case of these four demonstrations, typical of interferometric 
settings and instructive in figuring out how these circuits work. Suppose that all inputs to the 
circuit have the form (|2.11[) . with \ipo) and being orthogonal, but with a relative phase <j> 
between |?/>o) and Thus the inputs take the form 

gd|^0>®|0) + v^"e-^i)®|l). (2.24) 



We can apply all our results for the three circuits simply by making the replacement \tp\) — > 
e~ % ^\ip\). The phase <fi, which can be thought of as the phase to be estimated in an interferometer, 
is encoded in the entanglement of the upper wire of the circuit with the qubit on the lower wire. 
Furthermore, suppose that U swaps \ipo) and i.e., U\ipo) = and U\ipi) = \ipo). If we think 
of |V>o) and IV'i) as the standard states in a two-dimensional Hilbert space, then U is the Pauli X 
operator in this space, and P can be chosen to be the Pauli Z operator. 

In this situation it is clear that circuit (ii) employs a joint measurement to read out the phase 
information directly. In circuit (i), the role of the controlled-C/ is to disentangle the upper and 
lower wires, leaving the product state IV'o) ® {\/Qo\fy + y/0i e ~ lr ^ I !))> hi which the phase information 
is written on the qubit. The subsequent Hadamard transfers the phase information to amplitude 
information, which determines the outcome probabilities for the measurement of Z. Likewise, 
in circuit (iv), the state after the Hadamard and controlled- P is a product state, (^/%\ipo) + 
^/qle~ ^ ^\'^pl)) (g> (|0) + |l))/\/2, but in this case, the phase information is written onto the top line 
of the circuit. Circuits (i) and (iv) are thus both disentangling circuits, which write the phase 
information onto one of the quantum wires, where it can be read out by a measurement on that 
wire alone. For both circuits, the probability of outcome z is p z = |(1 + 2z^Jq§q\ cos 4>). 

Circuit (iii) is also disentangling, but works in a different way. With the current assumptions, 
the two outcomes of the measurement of Z are equally likely, i.e., r y = 1/2, and the state of the 
upper wire after the measurement of Z yields outcome y is y^lV'o) + y\^Qi^~ t ^\ 1 Pi}- The upper 
wire now carries the phase information, but one needs the outcome y to know how to extract phase 
information from the subsequent measurement of U. The conditional probability for the outcome x 
of the U measurement, given y, is p x \ y = ^(1 + 2xy^/qoqi cos <j>). If one does not use the outcome 
y, then the phase produces no interference, the unconditioned probability for x being p x = 1/2, 
but if one does use the phase information, as in the classical control of circuit (iii), the phase does 
produce interference, the unconditioned probability for the product, z = xy, being the same as for 
circuits (i), (ii), and (iv). 

In circuits (i) and (iv), disentanglement occurs as a consequence of coherent operations, but 
in circuit (iii), disentanglement is a consequence of the measurement of Z. This measurement- 



induced (or post-selected) disentanglement was originally called a quantum eraser [4l|, |42j, l43l ]: 
the modern understanding is that what is "erased" is the entanglement that destroys interference 
when one measures only the top wire or the bottom wire. A nice feature of the circuit (iii) is 
that one can regard the Hadamard and measurement of Z as occurring either before or after the 
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FIG. 7: Conventional qubit and equivalent modal interferometers, which use independent particles. The 
top circuit, a qubit interferometer, clearly acts independently on the input qubits, subjecting each to a tt/2 
pulse, a rotation that imposes the phase (f>, and a second tt/2 pulse. The state of the N qubits after each 
gate is given. Each qubit is measured in the standard basis; the probabilities to be found in states |0) and 
|1) are cos 2 (0/2) and sin 2 (0/2). The N outputs are processed to extract the difference 2m — no — n\. 
The middle circuit writes the same interferometer in terms of Schwinger operators, and the bottom circuit 
completes the translation to a modal description. The N photons, initially all in mode 0, are subjected to 
a 50/50 beamsplitter, a phase shifter, and a second 50/50 beamsplitter, after which the number of photons 
in each mode is counted, with the outputs differenced to produce 2m = no — n\. Quantum states are 
again tracked after each gate. The output 2m is binomially distributed; the mean and variance of m are 
(to) = t^N cos 4> and A 2 to = -|iVsin 2 0. The resulting phase sensitivity, 5<p — Am/\d{m)/d(j>\ = l/>/~N, is 
the quantum noise limit. The overall interferometric transformation is a rotation by —<j) about the x axis, 

i p piJyT /2 e -iJ z <t> p-iJy-rr /2 _ giJ x 



measurement of U. In the former case, they seem to prepare a state on the top wire that carries the 
phase information and in the latter case, they seem to "erase" the entanglement that destroyed the 
interference on the top wire, but there is no difference between these two cases, only a difference 
in perspective. 
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FIG. 8: Coherent-state version of a conventional modal interferometer. The input coherent state \a, 0) = 
D(a) (g> /|0,0) = e~l"l ^ 2 J2N=o( aN /v^-^N,®) is the result of applying the displacement operator D(a) 
to the vacuum state of mode 0. The displacement operator is the mathematical description of the action 
of a classical driving field on the mode; the action of an ideal laser operating far above threshhold is well 
approximated by a displacement operator, apart from a slow drift in phase. Since the interferometer itself 
is number-preserving, the different TV-photon sectors of the coherent state are processed independently, 
implying that the output probabilities are obtained by averaging the binomial distribution of Fig. [7J over 
the coherent state's Poisson distribution for total photon number. The resulting mean and variance of m 
are (m) = ^|a| 2 cos</> and A 2 m = ||a| 2 sin 2 0, which gives phase sensitivity 5(j) — 1 / |or| . We can think 
of the coherent-state interferometer as a classical interferometer contaminated by the shot noise of the 
coherent state's Poisson distribution. The classicality of the coherent states appears in the state description, 
tracked through the gates, as the fact that the phase <j> appears not as a relative phase between states in 
superposition, but as a phase change of coherent-state (classical) complex amplitudes. 



III. QUANTUM-NOISE-LIMITED INTERFEROMETERS 

The chief purpose of this paper is to examine and relate different methods for achieving 
Heisenberg-limited sensitivity, i.e., phase sensitivity that scales as 1/N with the number N of 
atoms or photons. Before taking up that task in Sec. IIVI however, we first look, in Fig. [TJ at 
conventional qubit and modal interferometers, which use independent particles. Qubit interferom- 
eters are typically atomic (Ramsey) interferometers, and modal interferometers employ two modes 
of the electromagnetic field. The sensitivity of a conventional interferometer scales as 1 /vN, the 
so-called quantum noise limit (or shot-noise limit) [l[. 

A couple of aspects of the interferometers in Fig. [7J deserve attention. In the modal version, 
the input Fock state \N) for mode can be replaced by a coherent state |a) = D(a)\0), where 
D(a) = e aa ~ a a is the modal displacement operator. It is worth recording that e~ lJ " nS D(a) (g> 
D(fi) e iJ-M = D(a') (8) D(j3'), where 

with M given by Eq. ()2.6|) . This implies that 

e - iJ -" e k/?} = |a',/3') . (3.2) 

A coherent-state interferometer is depicted and described in Fig. [8l It is a semi-realistic rendition 
of what is done in optical interferometry, since an ideal laser produces a very good approximation 
to a coherent state. 

It is clear from both the qubit and modal versions of a conventional interferometer that the 
output is insensitive to common-mode phase shifts. This is the primary advantage of measuring 
phase shifts in an interferometric setting. It is also clear, however, that the equatorial axes of 
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FIG. 9: (a) A measurement of Xn — \0)(N\ + \N)(0\ is equivalent to an effective Hadamard gate Hn = 
(Xn + Zn)/V% followed by a measurement of Zn — |0)(0| — |A)(A|. A subscript N on a two-dimensional 
unitary operator signifies that operator acting in the "qubit" modal subspace spanned by standard states |0) 
and | A). The final form, which replaces the measurement of Zn by photocounting, is valid if the input state 
lies in the "qubit" subspace spanned by |0) and \N); it is an example of extending a measurement outside a 
subspace of interest, (b) A measurement of X® N on N qubits is equivalent to measuring i N e~ lJmn on modes, 
and this in turn is equivalent to a 50/50 beamsplitter followed by a measurement of the parity of mode 1. 
Notice that X® N = i N e ~ iJ ^ is the qubit FLIP and modal SWAP operator: i N e~ iJ *' n \n , n x ) = K,ra ). 
We use this measurement in our discussion of NOON-state interferometers, where all that is really needed is a 
measurement of the observable |JV, 0)(0, N\ + 10, N) (N, 0| in the two-dimensional subspace spanned by \N, 0) 
and \0,N); i N e~ %J '°' K is a natural extension of this observable into the entire Hilbert space. Throughout 
this paper, we symbolize a parity measurement as in the last circuit — photocounting followed by classical 
extraction of parity — but the reader should remember that Fig.[6lc), with U = II, provides another, standard 
circuit for measuring parity, by mapping the parity to an ancillary qubit. 



the two 50/50 beamsplitters (or n/2 transitions) must be well denned relative to one another. A 
change in the equatorial axis of the second beamsplitter introduces a relative phase shift across the 
beamsplitter that is indistinguishable from the phase (p. In most applications, where one is only 
interested in changes in cj>, the actual requirement is that the relative phase shifts across the two 
beamsplitters be stable relative to one another over times somewhat longer than the time scale of 
changes in (f>. 

In different implementations, the stability requirements between the two beamsplitters assume 
different forms. In optical-frequency interferometers, where the modes typically correspond to 
different spatial paths or different polarizations, the phase shifts across the beamsplitters are de- 
termined by coatings on the beamsplitters and by placement or orientation of the beamsplitters; 
it is relatively easy to keep these stable. In atomic interferometers, the relative phase shifts across 
the 7r/2 pulses are chiefly determined by the times at which the pulses are applied; it is not difficult 
to keep the time separating the two pulses stable at microwave frequencies. 



IV. FOCK-STATE-SUPERPOSITION INTERFEROMETERS 



Entanglement among particles allows one to improve the sensitivity scaling for estimating a 
phase beyond that offered by the quantum noise limit. The optimal sensitivity scaling, known as 
the Heisenberg limit, is given by 1/N BQEi]. There are many different methods for achieving 
this limit, most of which involve using superpositions of two Fock states of one or two modes. In 
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FIG. 10: Four Hcisenberg-limited phase-estimation protocols that use Fock-state superpositions: (a) single- 
mode phase estimation; (b) single-mode, qubit-assisted phase estimation, (c) two-mode interferometer, called 
a NOON interferometer after the state before the phase shifter; (d) two-mode, qubit-assisted interferometer. 
In the qubit-assisted protocols, the qubit can be thought of as a second or third arm of an interferometer. 
Quantum states are tracked through each protocol. In protocols (b) and (d), inserting a gate e te / 2 e ~ lZ0 / 2 
between the two controls on the qubit wire shifts the fringe pattern by angle 9; in particular, an S gate shifts 
the fringes by n/2. In protocols (a) and (c), the output is either photons or X photons; mapping to +1 
and X to —1, the outcome becomes z = ±1, as in protocols (b) and (d). In all cases, the probabilities of 
the two outcomes are p z = ^(1 + zcosNcf). The mean and variance of z, (z) — cos N<f> and A 2 z = sm 2 N(j), 
give Heisenberg-limited phase sensitivity 5<j) = Az/\d(z) /d<fi\ — 1/X. The photocounting in the NOON 
interferometer (c) can be replaced by a measurement of Zn ■ Then the operations after the phase shifter arc 
equivalent to a measurement of Xn <8> Xn, which for the NOON-state input, is the same as a measurement 
of |X,0)(0, X| + |0,X)(X, 0|; this, in turn, as pointed out in Fig. [HJb), is equivalent on NOON-state inputs 
to a measurement of the modal SWAP operator i N e~ %J ^ . 
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FIG. 11: (a) Circuit equivalences that equate a controlled- SWAP between two modes to controlled-NOTs 
and controlled- SIGNs on N = 5 qubits, plus a Z N on the control qubit. (b) ^V-qubit cat-state interferometer, 
obtained by applying (a) to the two-mode, qubit-assisted interferometer in Fig. [TuTd). The interferometer 
is named after the cat state (|0)®^ A ' +1 ^ + |l)®^ Ar+1 ^ ) )/v / 2 at the input to the phase shifters; it has the same 
outcome probabilities and phase sensitivity as Fig. [TuTd) . Proposals for cat-state interferometers often also 
apply the phase shift to the control qubit; that does not happen here because of the special role of the 
control qubit, which is inherited from the two-mode, qubit-assisted circuit in Fig. UOf dV 



this section we consider four such methods and use circuit equivalences to show they are equivalent. 
In Sec. El we turn to corresponding methods that use superpositions of coherent states and find 
that they, too, saturate the Heisenberg limit, but with considerably better prospects for realization 
in the laboratory. 

Before getting started on Fock-state-superposition interferometers, we introduce two measure- 
ment rules in Fig. [9j The first, in Fig. EJa), is useful in our discussion of phase estimation using 
a single mode, and the second, in Fig. 0(b) , comes into play when we discuss genuine interferome- 
ters, which use two modes, such as NOON-state interferometers. In the first rule, we introduce the 
notation that a subscript N on a two-dimensional unitary operator denotes that operator acting 
in the "qubit" modal subspace spanned by standard states |0) and \N). The second rule is impor- 
tant for interferometers because it converts a measurement of the modal SWAP operator into a 
measurement that might be doable, i.e, a measurement of the parity of one mode behind a 50/50 
beamsplitter. 

The ability to convert to a parity measurement — and its importance — has been appreciated 



44 . |45| . Parity can be measured by photocounting followed 



m many papers 0, 0, n 0, m ei 

by classical extraction of the parity — this is how we represent a parity measurement throughout 
this paper — but it can also be measured using the protocol of Fig. [6tc) with U = II, which 
maps the parity to an ancillary qubit, and by other methods that rely on the interference of two 
modes 46j, |47], |48( . Indeed, these methods can be applied more generally to measuring the displaced 
parity operator, a measurement that comes up naturally in our discussion in Sec. W\ and that is 
the basis for direct measurements of the Wigner function of a mode [H, 2^, 3^, 40, 4a, 47, 48|, 49]. 
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FIG. 12: Circuit equivalences for disentangling circuits when the unitary U of Fig. [f^b) is X{f. Circuits 
are numbered as in Fig. [Hfb). Circuits (iii) and (iv) are converted to more useful forms by replacing the 
measurement of Xn by an effective Hadamard followed by photocounting, as in Fig. [Ha). In (iii), the 
classical controlled-EXCHANGE swaps the outcomes and N when the outcome of the Z measurement is 
y = +1. In (iv), Zjq is the operator that conjugates Xjq to its opposite. 



Figure [10] depicts circuit diagrams, labeled (a)-(d), for four phase-estimation protocols that 
employ superpositions of two Fock states. It is clear from the tracking of quantum states through 
these protocols that they are closely related and that all achieve Heisenberg-limited sensitivity by 
having a fringe pattern that oscillates N times as fast as the independent-particle fringe pattern 
of a conventional interferometer. We emphasize that these methods differ from conventional inter- 
ferometry in that they are not committed to placing beamsplitters immediately before and after 
the phase shifter, as in a conventional interferometer, but rather allow any procedure for preparing 
the input state to the phase shifter and any measurement procedure after the beamsplitter. 

Of the four protocols in Fig. [TOl the best known is interferometer (c), which inputs NOON 
states, {\N,0) dj0,i\A)A/2, to the differential phase shifter e~ iJ ^ and thus can be called a NOON 



, +J ( W N P| T 

interferometer Til . [r3 |. We remind the reader that the two-mode, qubit-assisted 

interferometer (d) is equivalent to an JV-qubit cat-state interferometer @, S, 0], as is shown in 

Fig. du 

The easiest equivalence to show is that of circuits (b) and (c). Insert a Z at the beginning of the 
qubit wire in (b). Push it to the end of the circuit through the measurement, thereby converting 
the controlled operations to control on |0). Replace the qubit with the two-dimensional modal 
subspace spanned by standard states |0) and \N). Insert an irrelevant global phase e lN ^^ 2 ; this 
changes the phase shift e~ m a< ^ on mode to a differential phase shift e~ lJz< ^ and completes the 
conversion of (b) to (c). 

The other equivalences, (b) to (a) and (d) to (c), are more interesting. It is easy to see that 
circuit (b) is a version of (a) with the phase information carried by an entangled state of the 
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FIG. 13: (a) Version of the circuit of Fig. [TOT a) obtained by applying the second form of circuit (hi) in 
Fig. [T2] to the circuit of Fig. nUTbh In this version, the superposition state input to the phase shifter is 
prepared by post-selection based on the outcome of the Z measurement on the qubit. The two possible 
outcomes occur with equal probability; the two input states, (|0) ± \N})/^/2, produce fringe patterns that 
are 7r out of phase. Averaging over outcomes destroys the fringe pattern, but the classical control circuitry 
restores a single fringe pattern. The initial circuitry and qubit measurement entangle and then disentangle 
the mode and the qubit; if the qubit measurement is regarded as occurring after the photocounting, this 
is the action of a quantum eraser, (b) Version of the circuit of Fig. [TPT a) obtained by applying the second 
form of circuit (iv) in Fig.[TJ]to the circuit of Fig. riuTb). Since Zn commutes with a^a, the controlled-Zjv 
passes unscathed through the phase shifter, becoming part of a coherent procedure for the superposition 
state input to the phase shifter. A final Hadamard gate is added to the qubit wire to return the qubit to its 
initial state |0). 



mode and the qubit; likewise, circuit (d) is a version of (c) with the phase information carried 
by entanglement of the two modes with the qubit. In the qubit-assisted protocols, (b) and (d), 
the post-phase-shift operations write all the phase information onto the qubit by disentangling it 
from the mode(s), whereas in (a) and (c), the phase information is on the mode(s). The key to 
showing the equivalences is thus to use the disentangling equivalences of Fig. [6^b) to move the 
phase information in (b) and (d) onto the mode(s) instead of onto the qubit. 

Figures [T2TI15I demonstrate the two equivalences; the demonstration that (d) is equivalent to (c) 
is very closely analogous to the demonstration that (b) is equivalent to (a) . Figures [12] and Q3] 
specialize the disentangling circuits of Fig. [UJ^b) to the situations in circuits (b) and (d), and 
Figs. [13] and [15] apply these to circuits (b) and (d) to convert them to versions of (a) and (c). The 
converted circuits have particular procedures, assisted by the qubit, for preparing an appropriate 
input state to the phase shifter. These preparation procedures are either post-selected, i.e., based 
on the outcome of a measurement on the qubit, or completely coherent. They require the ability 
to make A-photon states and to control two-dimensional Fock-state subspaces, both of which are 
very hard to do experimentally. Attempts to make NOON states beyond N = 2 up till now have 
had to be nondeterministic [121 ] . 

On the other hand, the superposition state that powers (a), (|0) + \N))/y/2, has been made 
deterministically for N = 1 to 5 in a qubit /cavity setting, where the qubit is a superconducting cir- 
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FIG. 14: Circuit equivalences for disentangling circuits when the unitary U of Fig. Efb) is e "TI ® 7 = 

^ F ®iV(_^®JV = x ®iv = iN e -u x -K = e - l j H 7 r / 2/ ^ ne4 j B 7 r /2_ The controlled- (n ® 7) is included to make 

t/ Hermitian as well as unitary. Circuits are numbered as in Fig. [6][b) . Circuits (iii) and (iv) are converted 
to more useful forms by replacing the measurement of i N e ~ lJ ^ by a 50/50 beamsplitter followed by photo- 
counting on mode 1 and classical extraction of parity, as in Fig.[9fb). The step to circuit (iv) is not generally 
valid, since it requires that the ±1 eigenspaces of i e" »* have the same dimension, which is only true 
when N is odd. We use circuit (iv), however, only when the modal inputs are confined to the subspace 
spanned by \N,0) and \0,N); in this case we can define P N by P N \N,0) = \N,0) and P N \0,N) = -\0,N). 
When N is odd, we can choose Pn = I ® n. 

cuit and the mode belongs to a microwave- frequency resonator [49( . The state-preparation protocol 
used in that experiment, due to Law and Eberly [50], can be used in principle to create arbitrary 
superpositions of resonator Fock states. It alternates rotations of the qubit with carefully timed 
Rabi couplings of the qubit to the resonator to deposit successive photons in the resonator mode. 
This protocol does not lend itself readily to the circuit notation of this paper, nor in its reliance 
on a ladder of two-level transitions, is it an efficient method for making large-iV superpositions 
(|0> + \N))/V2. 

Gerry and Campos [7[ have proposed a realization of the post-selected preparation procedure 
of Fig. [FoTa). The role of the qubit is taken over by two field modes that carry a single photon, 
and the controlled-e -4 ^ 71 ^ is achieved by a cross-Kerr coupling, which does a controlled swap of the 
coupled modes. 
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FIG. 15: (a) Version of the NOON circuit of Fig. [TOT c) obtained by applying the second form of circuit (iii) 
in Fig. [TH to the circuit of Fig. [TUf dV The controlled-(I <g) II) at the beginning of Fig. ITH i) is moved to 
the beginning of Fig. fTilf iii): it is pushed through the phase shifter and then through the controlled-e _l ' /! ' 7r , 
becoming a controlled- (I ® n) [because n ® Ie~ u ^ = (-Z)® N (-iY)® N = (-iY)® N Z® N = e^/^n], 
which can be omitted on mode l's vacuum input. In this version, the superposition state input to the phase 
shifter is prepared by post-selection based on the outcome of the Z measurement on the qubit. The two 
possible outcomes occur with equal probability; the two NOON inputs, (\N, 0) ± \0,N))/y/2, produce fringe 
patterns that are tt out of phase. Averaging over outcomes destroys the fringe pattern, but the classical 
control circuitry restores a single fringe pattern. The initial circuitry and qubit measurement entangle and 
then disentangle the modes and the qubit; if the qubit measurement is regarded as occurring after the 
photocounting, this is the action of a quantum eraser, (b) Version of the NOON circuit of Fig. flOTc) obtained 
by applying the second form of circuit (iv) in Fig. [14] to the circuit of Fig. fTOT d). The controlled- (II ® I) is 
handled as in (a). Since P/v commutes with J z , the controlled- Pm passes through the phase shifter, becoming 
part of a coherent procedure for the superposition state input to the phase shifter. A final Hadamard gate 
is added to the qubit wire to return the qubit to its initial state |0). The measurement procedure in (a) 
and (b) — 50/50 beamsplitter followed by photocounting and classical extraction of parity — doesn't look like 
the measurement procedure in Fig. flUf c) , but both procedures are equivalent to measuring the modal SWAP 
operator, i N 'e~ l,/a=7r as noted in Figs. [51 andfTUl 



V. COHERENT-STATE-SUPERPOSITION INTERFEROMETERS 



Just as for conventional interferometers, we can turn a Heisenberg- limited interferometer from 
one that uses superpositions of Fock states to a corresponding one that uses superpositions of 
coherent states. In this section, we take the single-mode and two-mode, qubit-assisted protocols 
of Fig. [TOTb) and (d) and convert them to run on coherent states. We then apply disentangling 
measurement equivalences to convert to protocols that prepare input states to the phase shifter 
that are superpositions of two coherent states. For a single mode, a superposition of two coherent 
states is often called a cat state, but we refer to it as a modal cat state to distinguish it from the 
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FIG. 16: Coherent-state analogue of the single-mode, qubit-assisted phase-estimation protocol in Fig. PTUT b). 
In constructing the analogue, one replaces the two controlled- A^s with controlled displacements that are 
inverses of one another. Here we make the first a controlled- D{— a) and the second a controlled-!? (a). We 
also introduce an S gate on the qubit to shift the fringe pattern by it/2 so that optimal sensitivity occurs 
at 4> — 0; this fringe shift turns out to fit naturally into further manipulations of the circuit. Quantum 
states are tracked through the circuit. An approximately equal sign (~), here and henceforth, indicates an 
expression that assumes small phase shifts, |a| 2 </> 2 -C \a\ 2 <fr -C 1. The approximate expressions emphasize 
the similarity to Fig. flOT b): exact expressions are given in the text. In the final state, we use the variable 
u = \a\ 2 <p — k/2 for brevity; it expresses the ir/2 phase shift produced by the S gate. For small phase 
shifts, the probability of outcome z is p z ~ i[l + z sin(|a| 2 0)], which gives Heisenberg-limited sensitivity 
54> = l/\a\. 

cat state of N qubits. 

Versions of our preparation procedures for superpositions of coherent states turn out to be 
practical, at least at microwave frequencies. Ultimately, there are two features that make coherent 
superpositions more practical than Fock-state superpositions. First, coherent states have their own 
phase, which we can play with in addition to the relative phase in superpositions; second, coherent 
states can be created and manipulated by coherent driving fields, represented by the displacement 
operator, thus providing an additional, crucial transformation that is not useful when dealing with 
Fock states. 



A. Single-mode coherent-state-superposition phase estimation 

Figure [16] gives the coherent-state analogue of the single-mode, qubit-assisted protocol of 
Fig. 110( b). The state after the controlled- D (a) is 

i=(|0)|0)+i J D(a)|-ae-^)|l)) . (5.1) 

It is easy to work out that 

D(a)\-ae-^) = D(a) e - iata(t, D(-a)\0) = e"^ 2 sin *| a (l - e"^)) , (5.2) 

from which one finds that 

(0\D(a)\-ae-^) = {a\e~^ a<t >\a) = e -l«l a (i-«-**) = e "*H 2 **4> e -2H* sin 2 (0/2) _ ^ 

In the small-angle approximation, |a| 2 2 <C \a\ 2 (j) ^ 1> t ne state ()5.2[) has negligible overlap with 
any Fock state other than the vacuum state and thus is given approximately by D(a)\— ae~ 1 ^) ~ 
e -A a \ ^|0). Written in a slightly different way, this becomes 

late-**) = e" iata *|Q) ~ e-^^a) . (5.4) 
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This is the approximation used in the state tracking of Fig. [T6j It says that a small phase shift of 
a coherent state simply places an appropriate phase in front of the coherent state; it expresses the 
fact that for the Fock states with appreciable amplitude in a coherent state, the phase shifts differ 
by roughly \a\cj) = \a\ 2 <p/\a\, which for small phase shifts is negligible compared to the average 
phase shift |a| 2 (/>. 

The probability for the Z measurement on the qubit to yield outcome z is p z = g(l + {%))■ The 
expectation value of Z in the final state is the same as that for X in the state (|5.ip before the final 
Hadamard. Thus this probability is given by 

1 
2 

-2|a| 2 sin 2 ((^/2)' 



Pz 



1 - zlm [{0\D(a)\-ae 



-i(j>\ 



1 + z sinl 



|q| 2 sin< 



-(l + zsin(|a| 2 </>)) 



(5.5) 



The final expression holds in the small-angle approximation and leads to Heisenberg-limited sensi- 
tivity for small phase shifts. 

A controlled displacement, like those in Fig. [TBI was introduced as a potential operation by 
Davidovich and co-workers 



13, 120], who called it a quantum switch, because a qubit coherently 



switches a classical driving field on a mode. These papers suggested that a quantum switch might 
be realized at microwave frequencies in a qubit/cavity set-up like that in Serge Haroche's Paris 
laboratory 23|, |2J]: a rubidium atom occupying one of two Rydberg levels, which constitute the 
qubit, passes through a superconducting microwave cavity; one level does not interact with the 
microwave field, but the other dispersively switches a cavity mode into resonance with a classical 
driving field, which excites the mode into a coherent state. 

Although a quantum switch has not been implemented in the laboratory, the same effect can 
be achieved by using a closely related procedure, depicted in Fig. [T71 which employs a sequence 
of displacements and controlled parities. A controlled-II is a particular case of a controlled phase 
shift, with the phase shift set to tt. Such controlled phase shifts have been implemented in the 
atom/cavity setting, with the dispersive coupling of one Rydberg level simply phase shifting the 
cavity mode. Such controlled phase shifts were proposed and analyzed in the atom/cavity setting 



in Ref. [181 ]. Initial experiments [23l . l51| had small phase shifts, but recent exp erimental work 



52j, [5J and vr [5J, |55|], 



in Haroche's group takes advantage of large controlled phase shifts, 7r/4 
produced in times much shorter than the cavity's damping time. The phase-estimation procedure of 
Fig. 1161 with the quantum switch replaced by the equivalent procedure of Fig. fTTl has been proposed 
and analyzed for implementation in the atom/cavity setting by Toscano and co-authors 56]. 

Our next task is to convert the qubit-assisted protocol of Fig. [16] to forms that have procedures 
for making modal cat states as input to the phase shifter. In this endeavor, Fig. [18] specializes 
the first three disentangling equivalences of Fig. [6](b) to the protocol of Fig. [TBI and Fig. [19] gives 
protocols in which the input state to the phase shifter is a modal cat state. 

The circuit in Fig. [T97 a) uses post-selection: it prepares one of two modal cat states, \ip y ) = 
(|0) + iy\a))/\/2, depending on which of the two equally likely outcomes y eventuates from the Z 
measurement on the qubit. The cat states are the relative states, \if)o) an d IV'i)) °f the extended 
discussion in Sec. Ill Bt they are normalized, but not orthogonal, having a purely imaginary inner 
product (ip-\ip + ) = i(0\a) = ie~^ I 2 . 

The protocol of Fig. [TWa) is precisely equivalent to that of Fig. [TBI since it is derived by using 
the disentangling equivalences of Fig. [18] to disentangle the qubit from the mode in Fig. [16] The 
state just before the photocounter, 



|^) = -L£>(-a/2)(|0} + 



iy\ae 



-i<t>\ 



(5.6) 
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FIG. 17: A quantum switch, in which a qubit coherently switches a classical driving field, is equivalent to a 
sequence of displacements and controlled parities, as one easily verifies using D(a/2)HD(— a/2) = D(a)H. A 
controlled-II is a particular case of controlled phase shift, with the phase shift set to it. The first controllcd- 
II can be omitted if the input state of the mode has even parity (e.g., the vacuum state). In our circuit 
diagrams, we often use the quantum switch because it is more compact, but the equivalent procedure is more 
realistic, at least at microwave frequencies, where the necessary controlled operations have been implemented 
in the laboratory. 

can be used to calculate the probability for outcome x of the parity measurement, p x \ y = |(1 + 
x(ip' y \n.\ij)y)). The expectation value of parity is 

<^|n|^> = i(<0|a> + (ae- < *|D(a)|-ae- i *)) - yIm(<0|D(a)|-ae-^>) . (5.7) 

We need 

(ae- l(t, \D{a)\-ae-^) = ^{ae'^aiX - e^)) = e -M 2 /2 e -4M 2 sin 2 «>/2) _ ^ 

Putting all this together, we get a conditional probability 

px{y = \{l + \xe-\ a \ 2 /\l +e -4|«| 2 sm 2 W2) ) +xysin( | a |2 sin0)e -2| Q Psin 2 W2 ) 

~ -(l + xe~ H2/2 + xysm(\a\ 2 (p)^ , (5.9) 

where the second expression assumes small phase shifts. As in the general discussion in Sec. Ill B[ 
averaging over y yields an unconditioned probability for x that has no ^-dependence, but the 
product variable, z = xy, reveals a fringe pattern, with the unconditioned probability p z given by 
Eq. 453]). 

The protocol of Fig. [TW b) uses a coherent procedure to prepare the cat state \tp+). It is easy to 
see that the probability for outcome z is given by Eq. (|5.9p with y = +1 and x = z. This probability 
is not the same as Eq. (j5.5l) . signaling that the circuit of Fig. [TW b) is not quite equivalent to Fig. [TBI 
though it becomes effectively so when \a\ is large, so that |0) and \a) are nearly orthogonal. 

Both circuits in Fig. [19] prepare modal cat states for input to the phase shifter, (a) by post- 
selection and (b) coherently. There is an independent coherent method for preparing a modal cat 
state, which does not require assistance from a qubit, relying instead on a self-Kerr coupling. This 
method, which goes back to the pioneering work of Yurke and Stoler [5t| |. follows mathematically 
from the simple identity e"™ 2 ^ 2 = e -WV(-i) n 7r/4 = e -"r/4[i + ^(-i)™]^ which implies that a 
self-Kerr coupling, with a tt/2 phase shift, makes a modal cat state: 

00 n -in 2 n/2 -i7r/4 

e -i(aUr./2 ]a) = e - H ' |n) = + <|-a» . (5.10) 



n=0 



The use of Kerr and higher-order nonlinearities to make cat states and more complicated superpo- 
sitions of coherent states has been elaborated extensively 58|, 5^]. 



Methods for preparing modal cat states by post-selection have been proposed using quantum 
switches 0, and using controlled phase shifts and controlled nonlinear interactions 0, [(3(], 
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FIG. 18: Circuit equivalences for disentangling circuits when the unitary U of Fig. [6jb) is Z3(a)II = 
D(a/2)UD(—a/2). The controlled-IT is included to make U Hermitian as well as unitary. Circuits are 
numbered as in Fig. |B£b), but circuit (iv) is omitted because it is more productive to guess a coherent 
preparation procedure to go with Fig. [16] than to determine an operator P that conjugates D(a)U to its 
opposite. Circuit (hi) is converted to a more useful form by replacing the measurement of D(a)II by a 
displacement followed by photocounting and classical extraction of parity. 



611 ] . and these have been applied to the cat-state interferometer [5] of Fig. [TWa). Dalvit and 
co-workers [(32| have investigated post-selected preparation of more elaborate superpositions of 
coherent states, such as superpositions of four coherent states, called compass states [63], for use 
in interferometers of this sort. The superpositions are created using quantum switches and Kerr 
interactions; the emphasis is on preparation procedures that might be implemented in an ion trap, 
where the qubit is two internal levels of an ion and the mode is a vibrational mode. In earlier work, 
Schneider and co-workers 64| reported methods for making such superpositions of coherent states 
in an ion trap. 



B. Two-mode coherent-state-superposition interferometers 



Having dealt with single-mode phase estimation, we turn now to the two-mode, qubit-assisted 
protocol of Fig. fTUT d). Figure [20] translates that Fock-state-superposition protocol into a corre- 
sponding coherent-state-superposition protocol. The translation replaces the preparation proce- 
dure of Fig. [I~07 d) with a procedure that uses two quantum switches to make the appropriate 
mode-qubit entangled state, as proposed in Refs. Iji, 20]. The rest of the circuit remains the same, 
except for the addition of an S gate on the qubit wire, which shifts the fringe pattern by n/2 
so that (f> = is an optimal operating point. This turns out to be a natural addition, just as in 
Sec. IV Al The final state provides the probability for outcome z, 
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l + zlm((a|e~ mta</ 
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■2|a| 2 sin 2 (0/2) 



(5.11) 



which gives Heisenb erg-limited sensitivity for small phase shifts. 
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FIG. 19: (a) Circuit obtained by applying the second form of circuit (iii) in Fig. [18] to the circuit of 
Fig. [HI The controlled-II at the beginning of Fig. [TST i) is moved to the beginning of Fig. [lSt iii), It is 
pushed through the phase shifter, S gate, and controlled- D(— a), converting the last to a controlled-D(a); 
it then can be omitted on the mode's vacuum input. The superposition state input to the phase shifter is 
prepared by post-selection based on the outcome y of the Z measurement on the qubit. The two possible 
outcomes occur with equal probability; the two input states, \tp±) = (\0) ± i\a)) / y/2 , are normalized, but not 
orthogonal modal cat states. The final state of the mode before photocounting is given in the small- angle 
approximation; exact expressions appear in the text. For small phase shifts, the conditional probability 
for parity outcome x is p x \ y — |(1 + xlj\) y ) ~ |[1 + xe~l Q I 2 + xy sm(\a\ 2 (f>)]. Averaging over y destroys 
the fringe pattern, but the classical control circuitry restores a fringe pattern for z = xy, with probability 
p z ~ |(1 + zsin(|a| 2 (/))], which gives Heisenberg-limited sensitivity S(f> = l/\at\. (b) Circuit that coherently 
prepares modal cat state as input to the phase shifter. Since D(a/2)UD(—a/2) swaps |0) and \a), 
we have D(a/2)HD(—a/2)\ip±) = ±i\ip T ). Thus we can disentangle the mode and qubit and prepare 
by inserting a controlled-II, preceded by D(—a/2) and succeeded by D(a/2), just before the qubit 
measurement in (a), which can then be omitted. Replacing the quantum switch by its equivalent from 
Fig. 1171 and adding terminal S and Hadamard gates to the qubit, thereby returning the qubit to its initial 
state |0), yields the circuit shown. We did not use circuit equivalences to obtain this coherent circuit, so 
it is not precisely equivalent to Fig. [TBI and to (a). For small phase shifts, the probability for outcome z is 
p z ~ |[1 + ze~' a ! / 2 + z sin(|a| 2 </>)], which reduces to the outcome probability in Fig. [TBI and in (a) when \a\ 
is large. 



The next task is to convert the qubit-assisted protocol of Fig. [20] to forms that have procedures 
for making the appropriate mode-entangled states, (\a, 0) ± i\0,a))/y/2, as input to the differen- 
tial phase shifter. The resulting circuits are given in Fig. [21] In CMC's research group at the 
University of New Mexico, we have taken to referring to these states as 0BB0 states (pronounced 
"oboe"), because we often write them (unnormalized) as \0,f3) + 1/3,0) [65J. We have analyzed the 
performance of 0BB0 states under losses and considered schemes for implementation, all of which 
is to be presented elsewhere (6ril ]. 

The circuit in Fig. [2TTa) uses post-selection: it prepares one of the two 0BB0 states, |\Py) = 
(\a, 0) +iy\0, a))/V2, depending on which of the two (equally likely) outcomes y emerges from the 
Z measurement on the qubit. The two 0BB0 states are normalized, but they are not orthogonal, 
their inner product, (\£_|\]/+) = i|(0|a)| 2 = ie~' Q ' , being purely imaginary. 
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FIG. 20: Coherent-state analogue of two-mode, qubit-assisted interferometer in Fig. [TOT d). To create the 
entangled state input to the differential phase shifter, the circuit uses two quantum switches, the first 
controlled on |1) and the second, via the surrounding X gates, effectively controlled on |0). An S gate on 
the qubit shifts the fringe pattern by n/2 so that optimal sensitivity occurs at <f> — 0. The entangled state 
input to the beamsplitter is a superposition of the states at the beamsplitter in Fig. llOf d). with coherent-state 
amplitudes in the superposition and with the extra i to shift the fringe pattern. Since the remainder of the 
circuit preserves photon number, the outcome probability for z can be obtained by averaging the probability 
in Fig. [TuT d) (with a ir/2 phase shift), |(1 — zsin A0), over the coherent-state Poisson distribution for N; 
the result is the distribution (|5. 1 1[) . which gives Heisenbcrg- limited sensitivity for small phase shifts. 



The circuit of Fig. [2lT a) is precisely equivalent to that of Fig. [20l since it can be derived from 
Fig. [2U] using the disentangling equivalences of Fig. [HI The state just before the photocounter, 



\#) = — ( 



ae-^ /2 /V2, -ae- i4,/2 /V2 ) + iy 



ae^ /2 /V2,ae^ /2 /V2)) , (5.12) 



provides the probability for outcome x of the parity measurement, p x \ y = |(1 + x($' \I U\ty' y )), 
via the expectation value of the parity of mode 1, 

IL\%) = (0\V2a) + lm((a/V2\e-^ a *\a/V2 } 2 ) 

= e"H 2 - y sm(\a\ 2 sin 0) e - 2 H 2 sin2 ^ . (5.13) 

Thus the conditional probability is 



Pxly = - + xe - |Q| -xysin(|a| 2 sin0) e - 2|Q| sin ^ /2) j ~ - (l + xe" |a| -xysm(\a\ 

2 ^ , - 2 (5,14) 

where the second expression is valid for small phase shifts. Although the unconditioned probability 

for x does not depend on <p, a fringe pattern appears in the product variable, z = xy, whose 

unconditioned probability is given by Eq. (|5.11|) . 

The corresponding protocol that prepares 0BB0 states coherently is depicted in Fig. I21lfb). The 

probability for outcome z is given by the conditional probability (|5.14|) with y = +1 and x = z. 

This probability is not the same as the final outcome probability in Fig. [20l telling us that this 

coherent protocol is not precisely equivalent to the circuits of Fig. [2D] and Fig. l2"TTa). It becomes 

equivalent when \a\ is large, so that \a, 0} and |0, a) are essentially orthogonal. 
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FIG. 21: (a) Circuit obtained by applying the second form of circuit (iii) in Fig. [14] to the circuit of 
Fig. HOI The controlled-II at the beginning of Fig. fTIT i) is shifted to the beginning of Fig. [T4T iii). It can be 
pushed through to just after the initial Hadamard on the qubit wire and can then be omitted on mode 0's 
vacuum input. The superposition state input to the phase shifter is prepared by post-selection based on 
the outcome y of the Z measurement on the qubit. The two equally likely outcomes produce input 0BB0 
states, \^±) = (|a,0) ±i|0, a)) /-\/2, which are normalized, but not orthogonal; they are superpositions, with 
coherent-state amplitudes, of NOON states, but with an additional i phase shift. Since the remainder of the 
circuit preserves photon number, we can calculate the probability for outcome z as a Poisson average of 
the NOON state probability of Fig. [ToT a) (with an additional ir/2 phase shift), p z = ^(1 — zsiniV0); the 
result is the probability p z of Eq. (|5.11|) . (b) Circuit that coherently prepares 0BB0 state as input to 
the differential phase shifter. Since D(a)H swaps |0) and \a), we can disentangle the mode and qubit and 
prepare |\&+) by inserting on both modes a controlled-II followed by a controlled-!? (a). Adding terminal S 
and Hadamard gates to the qubit, thereby returning it to its initial state |0), gives the circuit shown. Since 
we did not use circuit equivalences to obtain this coherent circuit, it is not precisely equivalent to Fig. [20l 
and to (a), though it becomes so when \a\ is large. 



The preparation procedure of Fig. [2"TYb). which employs quantum switches, is converted in 
Fig. [22Ta) to an equivalent protocol that uses controlled parities. This protocol is perhaps a prac- 
tical way to make an 0BB0 state at microwave frequencies because of the ability in an atom/cavity 
setting to do displacements and controlled phase shifts in times much shorter than the time scale 
for cavity damping. 

Figure [227 b) provides another, entirely independent coherent preparation procedure for 0BB0 
states, which was proposed by Gerry, Benmoussa, and Campos [9j. This procedure is a nonlinear 
interferometer, with a self-Kerr coupling in one arm. It generates an 0BB0 state by first making a 
modal cat state with the self-Kerr coupling and using the 50/50 beamsplitters of the interferometer 
to turn the cat state into an 0BB0 state. The nonlinear interferometer preserves photon number and 
thus can be converted to qubits using the equivalence of Fig. [2]^a) ; the resulting qubit interferometer, 
which relies on controlled-SIGN gates between all pairs of qubits, is depicted in Fig. [227c). The 
equivalent circuits (b) and (c) have been suggested and analyzed by Gerry and Campos [67J as a 
method for making a NOON-like state in a BEC. 
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FIG. 22: Coherent preparation procedures for OBBO states, (a) The preparation part of the circuit in 
Fig. I21f b), with quantum switches eliminated in favor of controlled parities by using the equivalence of 
Fig. [T7j (b) Kerr-gate protocol for making OBBO states. The circuit is a nonlinear interferometer, with 
50/50 beamsplitters before and after a Kerr self-coupling and a it phase shift on mode 1. (c) Translation 
to N = 3 qubits of the part of circuit (b) after the displacement of mode 0. The translation uses the 
circuit equivalence of Fig. EJa); the equivalence uses e~ lYlx l^ = HZ = XH. Since the part of (b) after the 
displacement of mode preserves photon number, if \a, 0) is replaced by \N, 0) at this point in the circuit, 
the output is the NOON-like state e _i7r/4 (|A^, 0) + i\0,N)/V2. Equivalently, (c) takes input \0)® N to the 
cat-like state e~"" r / 4 (|0)® iV + i\l)® N )/V2. On this input, the initial Zs can be omitted, and the final Xs 
only produce an overall phase, so they, too, can be omitted. 



VI. CONCLUSION 



After all this discussion of protocols for phase estimation, what can we say about the prospects 
for implementing Heisenb erg-limited phase estimation in the lab? Probably the best place to try 
is at microwave frequencies in a qubit/cavity setting, where coherent states can be manipulated 
relatively easily and where appropriate controlled operations are already available. One might 
think about performing one of the two-mode coherent-state protocols of Sec. IV Bl The qubit- 
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assisted protocol of Fig. [20l requires a controlled swap of the two modes, and the OBBO-preparing 
protocols of Figs. I2T1 and 1227a) require 50/50 beamsplitters and a measurement of parity. These 
are tall orders, except the parity measurement, which can be done, for example, by mapping to an 
ancillary qubit as in Fig. 0(c) . 

If one's primary goal, however, is Heisenberg-limited phase estimation — and not making 0BB0 
states — there is little reason to use two microwave modes. Almost certainly the most practical 
protocols are the single-mode coherent-state-superposition protocols of Sec. IV Al One can imagine 
implementing the qubit-assisted, single-mode phase-estimation circuit of Fig. [TBI with the quantum 
switches replaced by controlled parities, or perhaps the corresponding modal-cat-state protocols 
of Fig. [T9l with the quantum switches again removed in favor of controlled parities and the final 
parity measurement performed using an ancillary qubit. 

Despite the title of this paper, what we deal with here is not so much interferometry as general 
phase estimation in interferometric and other settings. Generally, we think about the phase shift 
as being applied to atoms, as in Ramsey interferometry, or to modes of the electromagnetic field, 
as in optical interferometry. 

In a conventional, quantum- noise-limited interferometer, there is a ir/2 pulse or beamsplitter 
immediately before and immediately after the application of the phase shift. Before the initial 
beamsplitter, there is some relatively simple state preparation, and after the final beamsplitter, 
there is a relatively simple measurement. Once we generalize to Heisenberg-limited phase esti- 
mation, however, we outgrow these beamsplitters. We must consider any preparation procedure 
before the phase shifter and any measurement procedure after it. In this more general setting, 
it is clearly better to use a representation that is not tied to the two beamsplitters and that can 
represent whatever goes on before and after the phase shifter. 

This is the reason that quantum circuits provide a far more useful way of representing gen- 
eral protocols for phase estimation and interferometry. Quantum circuits provide a platform- 
independent pictorial representation of all phase-estimation protocols, allowing us to switch easily 
from qubits to modes and to visualize — almost to experience, so tactile are the circuit diagrams — 
the equivalences between different sorts of interferometers and other phase-estimation protocols. 
The goal of this paper has been to develop tools for manipulating the quantum circuits for phase 
estimation efficiently. 
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